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Abstract 



The exact solution of the Cauchy problem for a generalized "linear" vectorial 
Fokker-Planck equation is found using the disentangling techniques of R. Feyn- 
Q\ ' man and algebraic (operational) methods. This approach may be considered as a 

generalization of the Masuo Suzuki's method for solving the 1-dimensional linear 
, Fokker-Planck equation. 

o 

00 

^ ; 1 Introduction 
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The Fokker-Planck equations (FPE), the one-dimensional FPE 
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-—[a(t,x)W] + — [D(t,x)W(t,x)}, t>0, xGR, 



dt dx ' dx 2 
and the "vectorial" FPE 



~dt 



-V • [a(t, x)w] + VV : {£>(*, x)w(t, x)} , i>0, xeR", (2) 



where a(t, x) = (ai(t, x), a 2 (t, x), . . . , a n (t, x)) T is the "drift vector", D(t, x) is a sym- 

d 2 D 

metric non- negative definite "diffusion" tensor field of II rank, and VV : D = dx _ d ^. 
(Einstein summation convention accepted), are widely used IE— [B5| as a tool in mod- 
elling various processes in many areas of the theoretical and mathematical physics, 
chemistry and biology as well as in the pure and applied mathematics and in engineer- 
ing: the nonequilibrium statistical mechanics (in particular in the theory of Brownian 
motion and similar phenomena: random walks, the fluctuations of the liquid surfaces, 
the local density fluctuations in fluids and solids, the fluctuations of currents, etc); 
the metrology ( Josephson voltage standards); the laser physics; the turbulence theory; 
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the cellular behaviour; the neurophysiology; the population genetics; the mathematical 
theory and applications of the stochastic processes, — to mention only a few of them. 

Because of its importance there have been many attempts to solve FPE exactly or 
approximately ( for a review see [4, 6 - 11, 14, 19]). Among the recent investigations 
on this problem noteworthy for us is the method of M. Suzuki |TJ| . 

In this paper we find the exact solution of following Cauchy problem: 



du 

— = ax(t)u(t, x) + a 2 (t) • Vu + a 3 (t)x ■ Vu + a A (t) : VVu, u(0, x) 



(3) 



where a^t) is a symmetric non-negative definite tensor function of second rank of the 
scalar parameter t. 

It is easy to see that the Eq.([3]) is connected with the "linear" vectorial FPE (|2]) 
with a linear in x "drift vector" a(£, x) = bi + 6 2 x and an independent of x diffusion 
tensor D. ( Here bi, b 2 and D are functions of t.) Therefore the Eq. (|3p is a slight 
generalization of the "linear" vectorial FPE (|2|) with t-dependent coefficients. 



In the paper [2C] the "isotropic" problems 

a\u{t, x) + a 2 ■ Vu + CI3X • Vu + a^Au, 



du 
~dt 



u(0,x) 



and 



du 

— = ai(t)u(t, x) + a 2 (£) • Vu + a 3 (£)x • Vu + aJt)Au, u(0, x) 
dt 



(4) 



(5) 



have been exactly solved ( here 04 and 0,4 (t) are arbitrary non- negative constant and 
function of t respectively). 

In the paper ]2T] we have found the exact solutions of the following Cauchy prob- 
lems: 

du 

— = aiu(t, x) + a 2 ■ Vu + o^x ■ Vu + : VVu, u(0, x) = 0(x) (6) 



and 

du 
~di 



dt 



ai(t)u(t,x) + a 2 (t) • Vw + a 3 (t)x- Vu + a 4 (t)a : VVu, m(O,x) = 0(x) (7) 



where S4 and a are symmetric non-negative definite tensors of second rank and 04 (t) 
is a scalar function; a 4 (t) > . ( It is obvious that the problem (JJ) is more general 
than the problem (0) : in (0), 64 (t) is arbitrary symmetric non-negative definite tensor 
function of second rank, while in (|7|) a 4 (t) has a special form: 0,4 (t) = a 4 (t)ci. ) 

Our method may be regarded as a combination of the disentangling techniques of 
R. Feynman |^2| with the operational methods developed in the functional analysis 
and in particular in the theory of pseudodifferential equations with partial derivatives 
23] — 1|27|| . As we have emphasized in [^0| and 0] this approach is an extension and 



generalization of the M. Suzuki's method |18| for solving the one-dimensional linear 
FPE®. 
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2 Exact Solution of the Cauchy Problem (|3|) 

In view of the t-dependence of the coefficients in the Eq. (JJ), formally we have for the 
solution of the initial value problem (j3|) an ordered exponential 

u(t,x) = (exp + ^ [a^s) + a 2 (s) ■ V + a 3 (s)x- V + a 4 (s) : VV] ds) 0(x), (8) 

where 



exp + / C(s)ds = T— exp / C(s)ds 



n=l 



1 + hm I dt i I dt 



tn-l 



dt n C{t 1 )C{t 2 )...C(t r , 



If we introduce the operators 

A(t) = a 2 (t) • V + a 3 (t)x • V and £?(£) = a 4 (t) : VV, 
we may write (§) in the form 



u(t, x) 



Jo ai ( s ) ds 



exp_ 



A(s) + B(s) dsU(x), 



(9) 



(10) 



11^ 



as the first term in the exponent commutes with all others. 

To proceed with the pseudodifferential operator in Eq. (|ll|) we shall use the theorem 
of M.Suzuki |§ : 

If 



then 



exp + 



A(s) + B(s) 



'A(t),B(tj\ = a(t,s)B(s), 
ds = (ex P+ J* i(s)ds) (exp + J* B(s)e~ So a M du ds 



In our case we have 

A(s),B(s')\ = [a 2 (s) • V + a 3 (s)x- V , a 4 (V) : VV 



= -2a 3 (s)a 4 (V) : VV = -2a 3 (s)B(s') 
Therefore from (|T2|) we obtain 

t 

'o 



exp_ 



i(s) + B(sj\ ds = (exp + J* A{s)ds\ (exp + ^' 5(s)e 2 /o a ^ du ds 



(13) 



(14) 



The linearity of the integral and the explicit form of A (see Eq. (|i~0"D ) permit to write 
the first factor in fll!]) in terms of usual, not ordered, operator valued exponent 

exp + J* A(s)ds = exp + jf [a 2 (s) ■ V + a 3 (s)x ■ V] ds = e^M^+^W*^. (15) 



For convenience we introduce the following notations: 



ai(t) = / ai(s)ds, a 2 (t) = / a 2 (s)ds, a 3 (t) = / a 3 (s)ds. 
J o Jo Jo 



Consequently (from now on '"" means 



at' 



a[(t) = ai(t), a 2 (t) = a 2 (t), a 3 (t) = a 3 (t), 

ai(0) = 0, a 2 (0) = 0, a 3 (0) = 0. 
Thus we obtain from the Eq. ( p]) 

u(t,x) =e °i(*) e [a2W+a 3 (*)x]-v ^exp + ^ < a 4 (s)e 2Q3(s) : VVds 
Finally using the formulae (see p(| and |2l[| ) 



x . 



exp+ (/ **(s) : VVds 
1 



exp 



where 



and 



^/det(47rf(t))^ 

dy = dyidj/2 • • • dy, 



e a 2 (t)-V+a 3 (t)x-V^/ x \ _ ^ ( xe «8(*) 



x-y) : (x-y) 



0(y)dy, 



f(t) = / m(s)ds 
Jo 



+ j\ 2 {s)e a ^ds^ =g( 



(16) 



(17) 



(19) 



(20) 



we find from the Eq. ([Tj|) the following expression for the exact solution of the Cauchy 
problem @ (\j/(s) = a 4 (s) exp[2a 3 (s)]) : 



it(£, x) 



3 ai(t) 



det(47rf(t))j£ 



exp 



-( z -y) • —4— • ( z -y) 



0(y)d?/, (21) 



where 



f(t)= f a 4 (s)e 2 " 3(s) ds 

JO 



is a symmetric non-negative definite second rank tensor function of t , dy — dyi . . . dy n 
and z is defined in (p0[). 

Substituting the expression (|2l|) in the Eq. @ we see immediately that the function 
u(t, x) is a solution of the problem (|3|), and, according to the Cauchy theorem, it is the 
only classical solution of this problem. 
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3 Concluding remarks 



• The exact solutions of the Cauchy problem (|J) is obtained using the algebraic 
method we have described. 

• When a 4 (t) is scalar: a 4 (t) = a 4 (t)l ( in this case a 4 : VV = a 4 A ) the 
"anisotropic" problem @ turns to the "isotropic" one, with the exact solution 
found in [EH]. It is easy to check that the solution ( plj ) turns to the solution 



obtained in pOf (there is an error in |j20 |: the sign before a 2 in the Eqs. (17) and 
(34) there, should be (+)). 

In the case a 4 (£) = a 4 (i)d the Cauchy problem (|3|) reduces to the problem (|7|) 
treated in [pjj. In this case the solution (pT|) turns to the solution obtained in ]2TJ . 



• For different choices of the coeficients aj and a 2 the Eq. (0) may be regarded 
also as a set of different diffusion equations. Therefore from the formula ( |2"T| ) 
we obtain the exact solutions of the Cauchy problems for this set of diffusion 
equations. 
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